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Abstract 

The non-perturbatively defined running quark mass introduced by the 
ALPHA collaboration is based on the PCAC relation between correlation func¬ 
tions derived from the Schrodinger functional (SF). In order to complete its 
definition it remains to specify a number of parameters, including the ratio 
of time to spatial extent, T/L, and the angle 9 which appears in the spatial 
boundary conditions for the quark fields. We investigate the running mass 
in perturbation theory and propose a choice of parameters which attains two 
desired properties: firstly the two-loop anomalous dimension df F is reasonably 
small. This is needed in order to ease matching with the non-perturbative 
computations and to achieve a precise determination of the renormalization 
group invariant quark mass. Secondly, to one-loop order of perturbation the¬ 
ory, cut-off effects in the step-scaling function are small in 0(a) improved 
lattice QCD. 


1 Introduction 


This paper is part of a project by the ALPHA collaboration to determine 
the A-parameter and the renormalization group invariant quark masses with 
controlled errors, using hadronic observables as experimental input []l|j2|. For 
numerical simulations of the lattice regularized theory, the basic difficulty con¬ 
sists in the large difference of length scales, ranging from the long distances 
typical for hadronic physics to short distances where perturbation theory can 
be applied with confidence. The proposed solution [Qj combines an intermedi¬ 
ate finite volume renormalization scheme with a finite size scaling technique J3[, 
which allows to step up the energy ladder recursively. 

The very definition of the renormalized running coupling and quark mass 
is of great importance for the method to be practical. The running param¬ 
eters should be relatively easy to compute by numerical simulation and they 
should not be affected by large cutoff effects. It is then possible to perform 
reliable continuum extrapolations and trace the non-perturbative evolution 
of the running parameters directly in the continuum limit, covering a wide 
range of energy scales. At high energies the evolution may be compared with 
perturbation theory, and, once perturbative evolution has set in, one may use 
perturbation theory to evolve to infinite energy and determine the renormaliza¬ 
tion group invariant parameters. For this last step to be feasible, the matching 
to perturbation theory should not require extremely high energies. This can 
be regarded as a further requirement to be met by a sensible definition of the 
running parameters. 

In this paper we carry out a perturbative investigation of a two-parameter 
family of running quark masses in the SF scheme [Q], which depends on the 
ratio p = T/L between time and spatial extent of the space-time manifold, as 
well as on the parameter 6 appearing in the spatial boundary conditions on the 
quark fields. In particular we determine the two-loop quark mass anomalous 
dimension df F which is needed for a precise determination of the renormaliza¬ 
tion group invariant quark mass ||]. It turns out that the parameters must 
be chosen with care for df F to be reasonably small. Taking the size of the 
one-loop cutoff effects in the step-scaling function as a further criterion sug¬ 
gests a specific choice of both parameters. This completes the definition of 
the running quark mass in the SF scheme which is expected to meet all of the 
above mentioned requirements. 

Many technical details of our perturbative calculation have appeared in 
ref. H and will not be repeated here. Parts of our results as well as of the cor¬ 
responding non-perturbative study ||] have already been published in ref. Q. 
This paper is organised as follows: In sect. 2 we recall some well-known facts 
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about the renormalization group, in particular the relation between running 
parameters and the renormalization group invariants. We then review the SF 
scheme and present our perturbative results for the renormalized quark mass 
in the continuum limit (sect. 3). Cutoff effects are discussed in sect. 4 and we 
end with a short summary. Finally, an appendix has been included to indicate 
the changes to appendix A of ref. § and appendix B of ref. || for our more 
general choice of parameters. 


2 Renormalization group 

In order to put the present work into its context we review some aspects of 
the renormalization group for QCD with N colors and A r f quark flavors and a 
diagonal quark mass matrix. 


2.1 Callan Symanzik equation and running parameters 

In the following it is assumed that the theory has been regularized, e.g. through 
the introduction of a space-time lattice, and that all renormalization condi¬ 
tions are independent of the quark masses. Any physical quantity P is a 
renormalization group invariant, i.e. as a function of the normalization mass 
(jl, the renormalized coupling g n and the renormalized quark masses m RjS , 
s = 1,..., Nf , it satisfies the Callan-Symanzik equation, 


f ,, o M 


d 


m H , s 


dm r 


P = 0. 


( 2 . 1 ) 


The renormalization group functions (3 and r can be derived from the rela¬ 
tion between the bare and renormalized parameters, using the independence 
of the former upon g. The precise formulae depend on the details of the 
regularization and will not be needed in the following. 

For small couplings (3 and r admit asymptotic expansions of the form 


OO 


P(g) s ~° —g 3 hg 2k , 

(2.2) 

k =0 


oo 

r{g) s ~° -g 2 ^2d k g 2k , 

k =0 

(2.3) 


with coefficients which are renormalization scheme dependent in general. In 
the minimal (MS) or modified minimal (MS) scheme of dimensional regular¬ 
ization b k and d k are known for k < 3 the first few being given by (with 
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Cf = (N 2 - l)/2 N), 


bo 

= {^AT - |Ar f }(4vr)- 2 , 

(2.4) 

d 0 

= 6C f (4tt)- 2 , 

(2.5) 

bi 

= {f N 2 - (f N- A- 1 )A f }(47r)- 4 , 

(2.6) 

di 

= C F {^f-N- fA" 1 - f AT f }(47r)- 4 . 

(2.7) 


While the minimal schemes of dimensional regularization are only defined in 
a perturbative framework, we emphasize that f3 and r are in general non- 
perturbatively defined functions. If these are given, running parameters at 
the momentum scale q are obtained by integrating the equations 




dm s 

1 ~di 


= r{g)m s , 


with the boundary conditions 

g(g) = g n , m s (g) = m R , s , s = l,...,N { . 


( 2 . 8 ) 


(2.9) 


The running parameters are related to the renormalization group invariant 
(RGI) quark masses M s and the A parameter, through (s = 1,... , Nf), 


A = q (b 0 g 2 ) bl/2b2 o exp j 

L 


1 


2&off 2 


x exp dx 


1 1 

+ 


bi_ 

(3(x) ' 6ox 3 bfa 

9 


M s = m s (2bog 2 ) d °/ 2b ° exp <j — J dx 


t(x) _ d 0 
(3(x) bpx 


( 2 . 10 ) 

( 2 . 11 ) 


One easily checks that the above expressions for A and M s are indeed solutions 
of the Callan-Symanzik equation ([2.1|). 

2.2 Finite renormalizations 

Any two mass independent renormalization schemes can be related by a finite 
parameter renormalization of the form 


I 1 = eg, 


c > 0, 

9r = 9r\J Xg(9n), 


m R, s = rn RjS X m (g R ), s = 1,... ,N { , 
3 


( 2 . 12 ) 

(2.13) 

(2.14) 












where we assumed the quark mass matrix to be diagonal in both schemes. The 
invariance of a physical observable P under such a change of variables implies 
the existence of a function P' such that 


P' {t 1 'in), 9 b.(9r), { m R, s (ffR) m R,s)}) = P((J,,g R ,{m RtB }), (2.15) 


and P' satisfies the Callan-Symanzik equation in the primed scheme with 
renormalization group functions /3' and t 1 . given by 


PisL) 

r'(g'n) 




9g'n 


dg b 


R ) ffR=9R(9k) 

d 


r(g R ) + /3(5r)w— In X m (g K ) 
@9 r 


9R — 9r(.9r) 


(2.16) 

(2.17) 


In perturbation theory, the finite renormalization constants X g and X m are 
expanded according to 


x{g* 


9r^' 

r^j 


1 + Y j X^g ] 


2k 

R ' 


(2.18) 


k =1 


Inserting into eqs. ( 2.11 ),( 2.17 ) one finds that bo, b\ and do are the same in 
both schemes (these are the “universal” coefficients), while all other coefficients 
are scheme dependent. In particular, the two-loop anomalous quark mass 
dimensions are related by 


d[ = di + 2 b 0 xW _ d 0 X 


(2.19) 


Renormalization group invariant parameters can also be formed in the primed 
scheme and one finds 


A' = A exp {xP/2b 0 }, (2.20) 

M's = M s , s = 1,... , Nf. (2.21) 


Therefore, given these parameters in some renormalization scheme, they are 
exactly known in any other scheme by just computing the one-loop relation 
between the renormalized coupling constants. Moreover, any physical ob¬ 
servable can be considered a function of the renormalization group invariant 
parameters, i.e. there exists a function P such that 

P( A, {M s }) = P{(x, g R , {m R , s }). (2.22) 


It thus appears natural to regard the non-perturbatively defined renormaliza¬ 
tion group invariants {M s } s =i.... ,jv f and A (in some renormalization scheme) 
as the fundamental parameters of QCD. 
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2.3 Determination of the RGI parameters 


In sect. 3 we will review the SF scheme, which is a non-perturbatively defined 
mass-independent renormalization scheme. The non-perturbative evolution of 
the running parameters in this scheme can be traced using numerical simula¬ 
tions and the functions g and fn s are then known up to some high energy scale 
q. One may look for the onset of perturbative evolution, and, once arrived in 
the perturbative regime, link the running parameters to the renormalization 


group invariants, using eqs. ( 2.1C ), ( 2.11 ) and the perturbative expansion of j3 
and r. In order to achieve the desired precision this requires the knowledge of 
the two-loop anomalous dimension df F and the coefficient b 2 F of the /3-function 
in the SF scheme Q. 

The easiest method to obtain df F consists in calculating the one-loop 
relation between the renormalized parameters in the SF scheme and some 
other scheme in which d\ is already known [cf. eq. ( 2.1 9|) 1 . For the purpose of 
our perturbative study the MS scheme is an appropriate reference scheme and 


we may thus use the result (2.7). 


3 The SF scheme 

We consider QCD on a finite (Euclidean) space-time manifold of size T x L 3 
with Schrodinger functional boundary conditions for the fields JlO] [E|. All 
dimensionful quantities are defined in units of L. In particular, for a given 
correlation function, the ratio p = T/L is assumed to be fixed so that there is 
(apart from quark masses) only a single scale, L, in the theory, which plays 
the role of the inverse normalization mass, L = 1/p j|]. 

In addition to dependence on the geometrical ratio p , correlation functions 
will also be functions of the parameter 9 appearing in the definition of the 
spatial boundary conditions on the quark fields (see e.g. eq.(4.8) in ref. [ 121). 
One is completely free to employ different choices of the parameters p and 9 
for the definitions of different physical quantities. This is a very convenient 
aspect of the SF framework which we will exploit; it can however lead to some 
slight abuse of notation which we hope is not confusing in our presentation. 


3.1 Renormalized coupling 

We start with the definition of the renormalized coupling constant in the SF 
scheme. This was introduced in ref. [10] for the pure gauge theory. The case 
of N = 3 colors, to which we restrict attention in this paper, was considered 


in detail in ref. [13], and we take over the background gauge field and choice 
of the parameter p = 1 specified there. 
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These choices were also made for the extension of the definition of the 


SF coupling to full QCD discussed in ref. [ 141. There the coupling was renor¬ 
malized at the physical value of the quark mass, thus leading to a quark 
mass dependent (3 function. Here we will deviate from this definition and 


adopt a mass independent renormalization scheme as advocated in ref. [ 121. 
Renormalizing the coupling at vanishing quark mass leads to the quark mass 
independent relation between the renormalized couplings 


9sf( L ) = 9M 5 (q)^g(9m( ( l)) 


(3.1) 


The SF coupling has been computed to one-loop order of perturbation theory 
in refs. 


1 


<Tg (1) = 26 0 ln(<?L) - —(ci, 0 + ciplVf) 


(3.2) 


where the coefficient cgi depends on the parameter 9. As a result of the 


detailed study in ref. [14] the particular choice 6 = 7r/5 was recommended for 
the definition of the SF coupling in QCD, in which case one obtains 


ci, 0 = 1.25563(4), c u = 0.039863(2). (3.3) 

In this paper we will also adopt this choice, but it will be clear where numerical 
results will change if in future simulations another choice of 9 is employed for 
g SF , and in fact it turns out that our general conclusions are not dependent 
on this. 


3.2 The role of the PCAC relation 

In QCD with JVf > 2 quark flavors the PCAC relation provides an attrac¬ 
tive starting point for the non-perturbative definition of a renormalized quark 
mass. Denoting the isospin non-singlet axial current and density by A“ and 
P a , respectively, the PCAC relation (for mass degenerate quarks), 


= 2 mP a 


(3.4) 


is a local relation between composite fields which is expected to hold when 
inserted in Euclidean correlation functions up to contact terms. 

The normalization of the axial current is conventionally fixed by requir¬ 
ing current algebra relations to assume their canonical form [jl5,16|. The non¬ 


linearity of these relations furthermore implies that the anomalous dimension 
of the axial current vanishes and so does the total anomalous dimension of 


the right hand side of eq. (3H). Therefore, a renormalized quark mass can be 
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defined through the PCAC relation by providing an independent renormaliza¬ 
tion condition for the axial density. Below the running quark mass in the SF 
scheme will be defined along these lines by using correlation functions derived 
from the QCD Schodinger functional. 


3.3 The renormalized axial density 

To define the renormalized axial density in the SF scheme we regularize the 
theory and choose the framework of 0(a) improved lattice QCD as discussed 
in ref. [|l^1. This choice is motivated by the corresponding non-perturbative 
studies Q, for which a perturbative investigation of the cutoff effects pro¬ 
vides complementary information (cf. sect. 4). We emphasize, however, that 
the results for the anomalous dimension are independent of this choice, and 
e.g. dimensional regularization would have been a practical alternative. In 
the following we use notations and conventions as in ref. (T^j without further 
notice. 

The renormalized 0(a) improved axial density has the form 

(■PnT = Zp(l + 6pam q )P a , P a = (3.5) 


where r a are the Pauli matrices acting in flavor space. If chosen appropriately 
the improvement coefficient bp cancels cutoff effects in on-shell correlation 
functions which are proportional to the subtracted bare quark mass m q = 
mo - m c . 

To define Zp we recall the definition of the bare correlation functions fp 
and /i, 

/p(®o) = -a 6 ^i(P a (x)C(y)7 5 ir a C(z)), (3.6) 

y,z 

h = -JJ 3(C'(u)75iP a C'(v)C(y)75^ a C(z))- (3.7) 

u,v,y,z 


The boundary source fields (, ( and (' are renormalized multiplicatively 
with a common renormalization constant Z£ [17,12], One may therefore define 
the renormalization constant of the axial density through the ratio 


Zp{go,L/a) = c 


Vfj 

MT/2Y 


(3.8) 


at vanishing quark mass m q = 0 and for vanishing boundary gauge fields. 
Here, the constant c is chosen such that Zp = 1 holds exactly at tree-level of 
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perturbation theory. Using the notation of refs. (6j,|j] with the modifications 
as indicated in the appendix we obtain 

c= ^L = \fN + 0{a 2 ). (3.9) 

v*o 


In general, c is a computable constant for a given lattice size and values of p 
and 6. 

The implicit dependence of Zp upon p and 6 will be discussed later. Here 
we emphasize that Zp is quark mass independent. Therefore, not only the (3 
function but also the anomalous dimension 


Tsf(.9r) — L 


dlnZp(g 0 ,L/a) 

dL 


90=3o(3r) 


(3.10) 


is quark mass independent. Here g R stands for g SF (L), and its (mass indepen¬ 
dent) relation to go is currently known to one-loop order |l4|] and to two-loop 
order in quenched QCD [18], 

In bare perturbation theory the renormalization constant Zp has an ex¬ 
pansion 


Zp(g 0 ,L/a ) = 1 + ^4 fc) (L/a)<7 0 2fc , (3.11) 

k =1 

(k) 

where in the limit a/L —y 0 the coefficients Zp ' are polynomials in ln(L/a) 
of degree k up to corrections of O (a/L). In particular the coefficient of the 
logarithmic divergence in Zp' 1 is given by the one-loop anomalous quark mass 
dimension do, and thus we parametrize Zp' 1 as 

4 1} = C F z p (d, p) - d 0 In (L/a) + O (a/L ). (3.12) 


Here we have made explicit the dependence of the cutoff independent term 
z p on 6 and p, which is inherited by any renormalized quantity involving the 
pseudoscalar density. The Feynman diagrams contributing to the correlation 
functions /p,/i (and hence to Zp) at one loop order have been discussed 
in detail in refs. H||. Here we use these results to extract Zp 1 for various 
different choices of the parameters p = T/L and 0. Some typical results for z p 
are collected in table 1. In section 4 we will also consider the 0,p dependence 
of the remaining cutoff terms in (3.12). 


3.4 Renormalized quark masses 


With the 0(a) improved bare axial current 

(M )m = + ca 4(5* + dp)P a , = ^7m75 


(3.13) 







p 

6 

z p (9,p) 

Zp |mom=0 

A z p 

di F /d 0 \ N=3 

2 

0.0 

-0.33879(1) 

-0.250 

-0.089 

0.9290 - 0.0441W 

2 

0.5 

0.07494(1) 

0.165 

-0.090 

-0.5880 + 0.0478Af 

1 

0.0 

-0.11953(2) 

-0.031 

-0.089 

0.1251 + 0.0046Af 

1 

0.5 

-0.09281(2) 

-0.003 

-0.090 

0.0271 + 0.0105JV f 


Table 1: The finite part of the one-loop renormalization constant of the axial density, the 
contribution of zero momentum gluon exchange, and the corresponding two-loop anomalous 
dimension. 


the renormalized current takes the form, 


(Ar)“ = Z A (1 + 6 A om q )(Ai)“. (3.14) 


To define the renormalized quark mass we first introduce a bare current quark 
mass through the PCAC relation between the unrenormalized fields, 


m = 


+ do)/ Al (xo) 

2/pOeo) 


x 0 =T /2 


(3.15) 


Here, the correlation function of the axial current is defined as /p in eq. ( p.6| ), 
but with the axial density replaced by the zero component of the (improved) 
axial current ( |3.13| ). 

The renormalized mass in the SF scheme is defined through the PCAC 
relation involving the renormalized 0(a) improved fields. This leads to the 
relation 


- rr\ (l + b A am q )Z A 

m SP (L ) = m——- = mZ A /Zp + 0(a). 

(1 + bpam q )Zp 


(3.16) 


In order to compute the two-loop anomalous dimension in this scheme we first 
relate rn SF (L) to the running mass in the MS scheme. To this end we start by 
combining the result for the axial current renormalization constant 1 19|, 161 


= -0.087344(1) x C F , 


(3.17) 


with the one-loop value of the ratio m q /m, which we obtained in the course 
of calculations done in ref. J5|. We then arrive at 

m SF (L) = m q {l + gl[d 0 ln(L/a) - (z p + 0.019458(1)) Cp\ + C%o)}- (3.18) 
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The corresponding relation for the renormalized MS mass has first been ob¬ 
tained in ref. | |19| and since then verified by many others, including one of the 
present author^. Then using the result 

m m (q) =m q {l+sg \~d 0 ln(aq) + 0.122282(1) xC F ] +0(^)}, (3.19) 


we obtain the one-loop coefficient 

= d 0 ln(qL) - (z p + 0.141740(2)) C F . (3.20) 


To finally obtain the two-loop anomalous dimension df F we may now use the 
known result for d\ in the MS scheme [cf. eq. (2.7)] and combine it with the 
one-loop coefficients (T2) and ( |3.2C| ) according to eq. Q. Proceeding in 
this way, we have avoided to expand eq. (|3.8| ) to order g 4 which would have 
required a two-loop computation. A few numerical values are given in table 1. 


3.5 Zero-momentum gluon exchange 

The results for Zp' 1 and d\ F with p = 2 show a strong dependence on the 
parameter 0. In the corresponding non-perturbative study a similar behavior 
is only seen at very small couplings Q. While this is not a problem in principle, 
it makes it more difficult to connect the perturbative regime to low energy 
physics along the lines of ref. [|j] • 

A closer look into the one-loop computation reveals that this strong in¬ 
dependence is almost entirely due to the exchange of gluons with zero spatial 
momentum [2C ]. These contributions are gauge invariant by themselves and 


it is not too difficult to compute them analytically. Setting 


w 


= Vsdp, 


(3.21) 


we find 

p 3 r 

z P {0, p) |mom=0 = — , -T 9 — \ 72 + 2Aw 2 - 5u> 4 - 96 cosh w 

192rc 4 cosh' w t 

+ (3u> 4 + 24) cosh 2 w — 12 w sinh 2rc|, (3.22) 

with the special case 


Zp(0,p)| mom =o = ~mP 3 - {3.23) 

1 Ref. |I^], although analytically correct, contains a small error in the quoted numerical 

results, which is caused by setting Poooi = 1.41 rather than Poooi = 1.310962... 
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One notices the overall factor p 3 which enhances the 0-dependence of z p for 
large p and results in the strong 0-dependence of df F for p = 2. If a subtracted 
constant is defined through 


AZp 


— Zp Zp | mom=0 5 


(3.24) 


the remaining 0-dependence is indeed very weak. We thus recommend to 
define the renormalized axial density with p = 1. Having fixed this parameter 
we furthermore choose 0 = 0.5 which leads to a conveniently small value of 
df F /do (cf. table 1). This completes the definition of the running parameters 
in the SF scheme. 


4 One-loop cutoff effects in the step scaling function 


An infinitesimal variation of Zp with the scale L defines the anomalous di¬ 
mension r SF (g R ) [cf. eq. ( 3. 10|) | . In the context of numerical simulations it is 
more convenient to consider finite variations of the scale, e.g. a change from 
L to sL, with a scale factor s. This leads to the definition of the step scaling 
function jlj 


£p (s,gl,a/L) = 


Z P (g 0 ,sL/a) 


Zp{go,L/a) 


90=3o(sr) 


(4.1) 


with continuum limit 


lim X P (s, g R , a/L) = op(s, g R ), (4.2) 

a —>0 

to be taken at fixed g R = g SF (L). 

We set the scale factor to s = 2 in the following. To one-loop order of 
perturbation theory we have 


£ P (2, gl a/L) = 1 + k{L/a)gl + O(gi) , (4.3) 

with 

k{L/a) = Zp ] (2L/a) - z£\l/o). (4.4) 


Using the notation of refs. ||||] for the renormalized correlation functions /p 
and /i, the one-loop coefficient on a finite lattice takes the form 


Z { p\L/a) 


t\ _ U\ _(i) ^2_ 

2to uq * 2fo 


U2 

Uo 


+ am. 


(i) 


£3 

2to 


^3 

u 0 


(4.5) 
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In order to see how fast the continuum limit 


k( oo) = —do ln(2), 


(4.6) 


is approached we define 


6^(L/a) = k(L/a)/k( oo) — 1. 


(4.7) 


This quantity contains all lattice artifacts at 0(g|). In the framework of 0(a) 
improved lattice QCD these are expected to decrease asymptotically with a 
rate proportional to a 2 / L 2 . In the present case, 0(a) improvement is achieved 
by setting Cgw =1- In particular, the boundary counterterms proportional to 
Ct and ct are not needed at this order of perturbation theory, owing, in the 
case of ct, to the identity 


t2 

2to 



(4.8) 


Some results for dp are tabulated in table 2. The numerical values have been 
obtained by inserting the exact expressions for the coefficients and 

for the given lattice size, and the coefficient am^ = —0.2025565(1) x 
Cf into eq. (|4.5|). 

Lattice artifacts appear to be reasonably small for all parameter choices 
considered. However, we observe that cutoff effects for p = 2 are generally 
larger than for p = 1. Furthermore the asymptotic O(a 2 ) decay seems to 
set in earlier for p = 2. It turns out that both effects are largely due to 
the zero spatial momentum gluon exchange contributions. Subtracting these 
contributions from dp we define 


A& = [k(L/a) - k(L/a)\ mom= o]/k(oo) - 1. (4.9) 


and list the numerical values in table 2. The subtraction terms are computed 
numerically for the given lattice size. In the special case 0 = 0 we also obtained 
a compact analytical formula, 

^/ O )L„„ 0 ,^(g-|S) C , (4.10, 

For p = 2 one clearly sees that the cutoff effects are dominated by the zero- 
momentum contributions. In the case 9 = 0 it is the explicit factor p in 
eq. ( |4.1C| ) which enhances the cutoff effects and also explains the early onset 
of the 0(a 2 ) behavior for p = 2. We conclude by noting that cutoff effects 
with the parameters p = 1 and 9 = 0.5 are indeed quite small, a fact that 
partially motivated this choice. 


12 






L/a 

O 

II 

Qb 

O 

II 

cb 

<1 

<$k 0=0.5 

A(5k 0 = 0.5 

P = 2 

4 

-0.3084 

-0.0025 

-0.2456 

-0.0702 

6 

-0.2292 

-0.0067 

-0.1499 

-0.0444 

8 

-0.1449 

0.0046 

-0.0893 

-0.0209 

10 

-0.0974 

0.0076 

-0.0584 

-0.0109 

12 

-0.0696 

0.0076 

-0.0412 

-0.0064 

14 

-0.0522 

0.0068 

-0.0307 

-0.0042 

16 

-0.0405 

0.0059 

-0.0239 

-0.0030 

P = 1 

4 

0.2040 

0.0650 

0.2136 

0.0140 

6 

-0.0121 

0.0126 

0.0208 

-0.0198 

8 

-0.0253 

0.0129 

-0.0026 

-0.0102 

10 

-0.0215 

0.0123 

-0.0062 

-0.0049 

12 

-0.0171 

0.0107 

-0.0064 

-0.0025 

14 

-0.0137 

0.0090 

-0.0058 

-0.0014 

16 

-0.0111 

0.0076 

-0.0052 

-0.0009 


Table 2: The one-loop cutoff effects <5k in the step scaling function, with and without zero 
momentum gluon contributions, for the same choice of parameters as in sect. 3. 


5 Summary 

We have carried out a perturbative investigation of a two-parameter family of 
running quark masses in the SF scheme. Its definition is based on the PCAC 
relation between correlation functions derived from the Schrodinger functional, 
together with an independent renormalization condition for the axial density. 
At one-loop order of perturbation theory and for asymmetric space-time vol¬ 
umes with T = 2L, many correlation functions show a strong dependence 
on the parameter 9. Its origin could be traced back to the contribution of 
gluon exchange with vanishing spatial momentum. Non-perturbatively this 
behavior is only matched at very short distances (J4|, making it more diffi¬ 
cult to apply the strategy of non-perturbative renormalization as outlined in 
refs. ]l|j2|. However, setting T = L completely eliminates this problem and 
also leads to a reasonably small two-loop anomalous dimension df F . We believe 
that this situation is generic and thus generally recommend the choice of sym¬ 
metric space-time volumes for the study of scale dependent renormalization 
constants in the SF scheme. 

In the present case we made the additional choice of 9 = 0.5 thus com- 
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pleting the definition of the running quark mass. One-loop cutoff effects in 
the step-scaling function are found to be reasonably small if 0(a) improved 
lattice QCD is used as a regularization. Finally, we mention that a corre¬ 
sponding non-perturbative study of the running quark mass is in progress Q, 
and preliminary results have been reported in ref. [|j. 

This work is part of the ALPHA collaboration research program. We 
thank Martin Liischer, Rainer Sommer and Hartmut Wittig for helpful dis¬ 
cussions and a critical reading of the manuscript. S. Sint acknowledges support 
by the U.S. Department of Energy (contracts DE-FG05-85ER250000 and DE- 
FG05-96ER40979). 
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Appendix A 


The expansions of the functions ti,Ui,Vi,Wi,yi i / 1 up to corrections of 
order ( a/L ) 2 for general p = T/L are as given in appendix B of ref. [Bl] and in 
appendix A of ref. ||] except that si, co are now defined by 

si = sinh(\/3#p), (A.l) 

co = cos h(V3 9p), (A.2) 


and the expressions for i = 3 now read 
t3 = 


2 A si L 2 N 

+ 


V3 = 


\/3#co 3 a co 2 

r 19A0si ApN9 2 (l - 2si 2 ) 

L 3\/3co 3 


3 co 4 


1- 

)V 


N si L 


U3 = 


V3 = 


r NO si 2 pN9 2 , 2,10 


+ 


\/3#co 2 a i6\/3co 2 

Asi(co — 2) L jA0 si(co — 2) 

y/39 co 3 a 

2 pN9 2 


3 co 4 


^ 6\/3co 3 

(co 3 — 4co 2 — 2co + 6) j-^-, 


w 3 = 


2 N L | 3 NO 2 + 4pA0 3 si j a 


+ 


CO a l co \/?>co 2 ' L* 


Asi(co + 2)L rA0si(co + 2) 


+ 


3\/39co 3 a 
2pN9 2 


^ 18\/3co 3 


9 co 4 


(co 3 + 4co 2 — 2co — 6) > — 


L 


(A.3) 

(A.4) 


(A.5) 

(A.6) 


(A.7) 
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